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^ ' ABSTRACT 

m ; 

^5 Context. Dust clouds influence the atmospheric structure of brown dwarfs, and they affect the heat transfer and change the gas-phase chemistry. 
However, the physics of their formation and evolution is not well understood. The dust composition can be predicted from thermodynamical 
equilibrium or time-dependent chemistry that takes into account seed particle formation, grain growth, evaporation and drift. 
' Aims. In this letter, we predict dust signatures and propose a potential observational test of the physics of dust formation in brown dwarf 
Or atmosphere based on the spectral features of the different solid components predicted by dust formation theory. 

q . Methods. A momentum method for the formation of dirty dust grains (nucleation, growth, evaporation, drift) is used in application to a static 
'_h ' brown dwarf atmosphere structure to compute the dust grain properties, in particular the heterogeneous grain composition and the grain size. 
^ . Effective medium and Mie theory are used to compute the extinction of these spherical grains. 

■ Results. Dust formation results in grains whose composition differs from that of grains formed at equilibrium. Our kinetic model predicts 
^ ' that solid amorphous Si02[s] (silica) is one of the most abundant solid component followed by amorphous Mg2Si04[s] and MgSi03[s], while 
Si02[s] is absent in equilibrium models because it is a metastable solid. Solid amorphous Si02[s] possesses a strong broad absorption feature 
centered at 8.7 /im, while amorphous Mg2Si04[s]/MgSi03[s] absorb at 9.7 jim beside other absorption features at longer wavelength. Those 
5— i features at A < 15yum are detectable in absorption if grains are small (radius < 0.2 fim) in the upper atmosphere as suggested by our model. 

Conclusions. We suggest that the detection of a feature at 8.7 jim in deep infrared spectra could provide evidence for non-equilibrium dust 
formation that yields grains composed of metastable solids in brown dwarf atmospheres. This feature will shift towards 10/im and broaden if 
silicates (e.g. fosterite) are much more abundant. 
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1. Introduction 



Atmospheric models are essential in the interpretation of com- 
plex brown dwarfs spectra, which are dominated by strong 
molecular absorption lines. The physics and chemistry of sub- 
stellar objects (brown dwarf and giant planets) are more com- 
plicated than anticipated because of non-equilibrium processes 
(dust formation, convective mixing) and their non-linear feed- 
backs on the radiative transfer. Consequently, parametrisa- 
tions of processes like dust formation and convection were ap- 
plied in atmospheric models so far. Current model atmospheres 
provide extensive solutions of the radiative tra nsfer prob- 
lem for the gas phase in hydrostatic equilibrium JTsuii et all 
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2001). However, the presence of dust is at present only in- 
ferred from classical models, where the dust composition 
is derived from gas/solid phase equilibrium considerations 
and time -scale arguments, and their comparison to observed 
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2003). Such equilibrium chemistry results in grains with a ho- 
mogeneous composition. 

Considering the non-equilibriu m character of phase tran- 
sitio ns (supersaturation» 1) IWoitke & Hellind J20031 12004) 
and lHelling & Woitkel J20061) proposed a theoretical approach 
to consistently model the formation of dust grains by nucle- 
ation (seed formation), growth, evaporation and drift (gravita- 
tional settling). In contrast to the phase-equilibrium calcula- 



2 



Ch. Helling et aL: Detectability of dirty dust grains in brown dwarf atmospheres 



tion, Woitke & Helling argue that grains in the oxygen-rich 
environment of a brown dwarf are heterogeneous in chemi- 
cal composition and in size. Moreover, precipitating into the 
denser inner atmospheres, these particles can reach sizes of 
several 100 fim. In Helling & Woitke (2006), the inferred dust 
composition differs markedly from that predicted by equilib- 
rium chemistry. These differences manifest in the intrinsic ab- 
sorption signatures of the solids in the mid-infrared range. 

This letter suggests a potential observational test for the 
presence of dust. The dust may be present in form of cloud- 
like structures. We will focus on the main characteristics of the 
theory, where the predictions in terms of dust spectral features 
can be tested with current and future instruments. In Sect. 13 we 
report on recent progress in modeling the formation of chemi- 
cally non-homogeneous cloud layers made of dirty (i.e. a vari- 
ety of compounds) grains in brown dwarf atmospheres. A de- 
tailed Mie theory treatment combined with effective medium 
theory allows us to suggest possible spectral features of such 
a dust layer for both approaches. The result of the modeling is 
described and discussed in Sect. [3] 

2. Approach 

We model nucleation, heterogeneous growth, evaporation, and 
drift (gravitational settling) of dirty dust p articles in a quasi- 
static atmosphere by the moment method jHelhng & Woitkel 
2006). Our prescription considers the formation of compact 
spherical grains in an oxygen-rich gas by the initial nucle- 
ation of Ti0 2 seed particles followed by the growth of a 
dirty mantle made of various solids. We consider amorphous 
Ti0 2 [s], Si0 2 [s], MgO[s], MgSi0 3 [s], Mg 2 Si0 4 [s], Al 2 3 [s], 
and metallic iron Fe[s], assuming that silica and silicates have 
approximately the same sticking probability. The moment and 
elemental conservation equations for Ti, Si, Mg, Fe, Al, and 
O are evaluated on top of a prescribed static model atmo- 
sphere structure (Allard et al. ') without further iteration on 
the temperature profile. The model calculates the amount of 
condensates, the mean grain size (a), and the volume frac- 
tions V s of each material as function of height z in the brown 
dwarf atmosphere for given model structures. At each depth 
in the atmosphere, the volumes fractions are used to com- 
pute mean dust optical constant s using the Maxwe ll-Garnett 
effective medium theory dBohren & Huffman! 1 1983k which is 
a valid method because the grains generally have one main 
volume component. The main dust component is the matrix 
in which the spherical inclusions made of the minor com- 
ponents are embedded. Dust extinction coefficients are sub- 
sequently computed using a Mie theory code for compact 
spherical grains. The different components are assumed in 
an amorphous state, consistent with laboratory experiment s 
of heterogeneous dust formation (e.g., Rie tmeiier et alJ fl999>. 
The optical constants for the amorphous s olids are o btaine d 
from la boratory measurements (TiC hfsl: iRibarskvl dl985h. 
SiQ2r f i: Eenning & Mutschkd i ll 9971 MgQ : jHofmeister et alJ 
12003). MgSiO^sl/MgoSiOJslT ljager et alJ d2003b. A l 2 (Ms |: 
iBegemann etafl i 19971) . and Fefsl: IOrdaletal.1 dl985ln . 

1 adopted from ftp.ens-lyon.fr/pub/users/CRAL/fallard/ 



Absorption cross sections and optical depth are then derived us- 
ing the dust density profiles assuming a grain size distribution 
f(a,z) = (a(z))6(a - (a(z))). For comparison, similar compu- 
tations were performed assuming that no Si0 2 [s] can form by 
setting the formation rate to zero. We also estimate the possible 
absorption depth in observed spectra. The contrast between the 
depth of the features and the continuum is an important param- 
eter for predicting the detectability of the dust features. 

3. Results 

We first show the grain size and dust composition profile of a 
AMES (cond) model with T eff = 1500 K, log g= 5.0 ( : ) and 
solar element composition. This can be considered a typical 
late L-dwarf atmosphere model. The outputs are then used to 
compute the dust extinction coefficient in the mid-infrared. 

3.1. Global grain size distribution and material 
composition 

The results depicted in Fig.Hlsho w the same glo bal dust cloud 
structure as presented in Woitke & Helling l J2004I) . The forma- 
tion of the cloud is governed by the hierarchical dominance of 
nucleation (upper most layers), growth & drift (intermediate 
layers), and evaporation & drift (deepest layers). The upper- 
most cloud layers are predominantly filled with small grains of 
a mean size of 10~ 2 /im which grow on their way into the atmo- 
sphere to 200 . . .300jum (l.h.s. Fig.[0. Eventually they enter 
even hotter atmospheric layers where they are no longer ther- 
mally stable. Hence, the grains shrink in size and finally dis- 
solve into the surrounding hot and convective gas. This picture 
reflects the stationary character of the grain component forming 
the brown dwarf's dust cloud where dusty material constantly 
fall inward and fresh, dust-free material is mixed upward. 



Table 1. Dust grain's mean size and material composition. 



local 


mean 


grain 


material 


temperature 


size 




composition 


T[K] 


{a) [fim] 


[volume frac] 


< 1000 


io- 2 


. . . io- 1 


25 - 35% Si0 2 [s] 
20 - 25% Mg 2 Si0 4 [s] 
10 - 25% MgSi0 3 [s] 
15% Fe[s], 10% MgO[s] 


1000 . . . 1700 


lo- 1 


...30 


25 - 30% Mg 2 Si0 4 [s] 

20 - 30% MgSi0 3 [s] & Si0 2 [s] 

15% Fe[s], 10% MgO[s] 


1700 .. . 2000 


30 .. 


. 200 


strongly changing composition 


2000 . . . 2300 


200 . 


.. 300 


80%Fe[s], ~ 18% A1 2 3 



The r.h.s. of Fig.^demonstrates the chemical composition 
of the dust grains. The grains are not of a single, homogeneous 
material composition and their material composition changes 
on their descent into the atmosphere. Since gas and dust ther- 
malise faster than the dust growth/evaporation processes take 
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Fig. 1. Left: Mean grain size stratification for a plan-parallel AMES (cond) late L-type brown dwarf model atmosphere of 
T e ff = 1500 K, log g= 5.0 of original solar abundance. Right: Chemical material composition of the grains in volume fractions of 
the solids for the same model (V^f- = £ s V s , s - contributing solid materials). Note the upside-down scale of the local temperature 
T[100K]. 



place (Woitke &Hellingll2003h . the chemical composition is 
determined by the local temperature and the reaction kinetics. 
Therefore, the upper atmosphere is populated by small, silica- 
and silicate-like grains until these materials become thermally 
unstable. Therefore in the lower and hotter part of the atmo- 
sphere big grains appear which are merely made of iron and 
some impurities of A^C^s] and TiC^ts]. The main dust com- 
ponent in the upper part of the dust layer is amorphous SiC<2[s] 
followed by Mg 2 Si0 4 [s]/MgSi0 3 [s] (see TableQ). The predic- 
tion contrasts with those from equilibrium dust formation mod- 
els where metastable species such as SiC>2[s] cannot exist (e.g 
Lodders & Fegley 2005). 

3.2. Spectral dust cloud features 

Figure |3 depicts the resulting spectral features from 6 /mi . . . 
15 fim for a brown dwarf dust cloud layer with the mean grain 
size and chemical dust material composition shown in Fig.[Q 

For grain composition predicted by heterogeneous dust 
formation, the only dust features with an appreciable con- 
trast are those of SiChts] centered at 8.7 fim and of 
Mg2SiC»4[s]/MgSi03[s] at 9.7 fim, with weaker absorption fea- 
tures around 20 fim and 32 fim (not shown). The features 
are broad (~ 1-1.5 fim) and lack substructures because the 
grains are amorphous. The abundance of metallic iron is high 
(~ 15%) but metals absorb photons continuously and do not 
show spectral features. The abundance of Al203[s] and TiC^ts] 
in grains are too small to significantly affect the overall ex- 
tinction coefficient. In models where SiC>2[s] is disregarded 
(dashed line, l.h.s. Fig. [2), the dust extinction is dominated by 
Mg2SiC»4[s]/MgSi03[s] at 9.7 fim. The mean grain radius re- 
mains smaller than 0.2 fim in the Td ust < 1 region in the model. 
Therefore, the grain absorption features remain. In contrast, if 
grains were a large as 10 fim, the resonance features would dis- 
appear jMin et all 120041) . The region where Td us t < 1 in the 



wavelength range 5-25 fim has pressure below 10 2 bar in the 
model atmosphere considered. 

In phase-equilibrium, the condensates are formed locally 
and retain the energetically most favorable composition at 
the temp erature and pressure of the given depth in the atmo- 
sphere (Lodders & Feglev, 2005 . in pressh . Further condensa- 
tion is not possible because the grains remain in relative thin 
discrete cloud layers and, hence, can not move into regions 
with still favorable growth conditions. Individual cloud layers 
are composed of different solid species depending on the tem- 
perature of the layer. Above the silicate clouds, the gas phase is 
strongly depleted and thus the absence of gas-phase species in 
L and T dwarf atmosphere has been used as an evidence for the 
presence of clouds. The prediction of equilibrium theory is the 
absence of solid SiC<2, hence equilibrium-dust is predominately 
composed of silicates. We modeled the absence of SiC<2[s] and 
the spectral result is shown in Fig.|2](l.h.s.). Beside moving the 
peak of the absorption to 9.7 fim, the absence of amorphous 
Si02[ S ] decreases substantially the absorption contrast. 

3.3. Discussion of detectability 

We now discuss the detectability of infrared dust features in 
brown dwarf atmospheres. From Fig.|2](r.h.s.), we can see that 
the absorption over continuum contrast varies from 10 . . . 20% 
with the maximum obtained for cold brown dwarfs if SiC>2[s] is 
present. This contrast is achieved because the grains are small 
(< 0.2 fim). In the absence of Si0 2 [s], Mg 2 Si04[s]/MgSi0 3 [s] 
will absorb weakly at 9.7 fim. In the case of large grains, no 
feature can be detected and it will not be possible to determine 
the composition of these grains. The presence or absence of 
SiOjTs] can be tested observationally. Recently, ICushing et alJ 
(2006) presented Spitzer IRS 2 results showing L-dwarf spectra 
with a 1 - 1.5 fim broad absorption centered at A » 9 fim. 



R = 90forA = 5.3/xm . . . 15.3 fim 
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Fig. 2. Absorption transmission spectrum (= F v /F cont with 
F cont the local continuum flux with respect to the feature) for 
the model in Fig.[0with (solid line) and without (dashed line) 
formation of SiC>2[s]. The absorption feature of SiChts] cen- 
tered at ~8.7 /an. 

However, we have treated an idealized situation where dust 
grains are the main sources of opacity. In reality, the feature 
will appear on the top of gas absorption lines. The gas phase 
absorption proceeds by lines whose width is determined by the 
pressure at the height of the absorption layer. Typical molecular 
band widths are ~ 0. 1 fim. 

Another limitation of our treatment is the assumption of 
spherical grains which may introduce spurious effects during 
the opacity calculation. Grains are most likely non-spherical 
and their absorption coefficients shoul d be computed with, for 
example, the hollow sphere model (Minetal., 2005) which 
successfully models asphericity by applying a distribution of 
hollow spheres. However, the use of hollow spheres will in- 
troduce minor changes in the shape of the absorption feature 
and will not change significantly the wavelength-position of the 
dust features and hence the conclusion of this paper. Finally, 
the heterogeneous formation model does not compute the ex- 
act grain size distribution but only its moments from which the 
mean size can be computed. Future work is needed to calculate 
better analytical representations of the size distribution func- 
tion from the dust moments. 

4. Conclusions 

We have studied the dust spectral signatures predicted by the 
Helling & Woitke (2006) models which make specific observa- 
tional predictions in terms of grain composition and size at dif- 
ferent altitudes. Heterogeneous dust formation in brown dwarf 
atmosphere results in dusty grains composed mainly of amor- 
phous SiChts] and amorphous Mg2Si04[s]/MgSi03[s] with 
impurities of Fe and Al oxides, TiC>2 and MgO. The dust forms 
one continuous layer in which the grains gradually change 
composition and size. The grains are sufficiently small in the 
upper atmosphere to produce solid res onance feature s in th e 
infrared, a conclusion reached also by ICushing et alJ J2006I) . 



In contrast, grain composition derived from chemical equilib- 
rium models is dominated by Mg2Si04[s] with the noticeable 
absence of amorphous SiC>2[s], which is a metastable species. 
The difference in dust composition allows us to make predic- 
tions in dust spectral features in the infrared for the two types 
of grains. Amorphous SiC>2[s] moves the peak of the absorption 
to 8.7/mi while a nearly pure Mg2Si04[s]/MgSi03[s] absorbs 
at 9.7 /j.m with a much smaller contrast. Therefore, infrared fea- 
tures of dust grains could be used to directly detect the presence 
of dust in brown dwarf atmospheres. The detection of dust ab- 
sorption features supports the prediction that grains have small 
size < 0.5 yum and is present in the Td ust < 1 upper atmosphere. 
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